Optimizing thermal transport in the Falicov-Kimball model: binary-alloy picture 



(N 
O 

o 

> 
O 

OO 



i 

CO 

-4-5 
c3 



I 

C 

O 

o 



> 

in 
oo 
m 



(N 
O 

c3 



o 
o 



X 



J. K. Freericks 1 and V. Zlatic 2 
1 Department of Physics, Georgetown University, Washington, DC 20057-0995, U.S.A. 
2 Institute of Physics, Zagreb, Croatia 
(February 1, 2008) 

We analyze the thermal transport properties of the Falicov-Kimball model concentrating on locating 
regions of parameter space where the thermoelectric figure-of-merit ZT is large. We focus on high 
temperature for power generation applications and low temperature for cooling applications. We 
constrain the static particles (ions) to have a fixed concentration, and vary the conduction electron 
concentration as in the binary-alloy picture of the Falicov-Kimball model. We find a large region of 
parameter space with ZT > 1 at high temperature and we find a small region of parameter space 
with ZT > 1 at low temperature for correlated systems, but we believe inclusion of the lattice 
thermal conductivity will greatly reduce the low-temperature figure-of-merit. 
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I. INTRODUCTION 

There has been a recent resurgence of interest in solid- 
state devices for thermoelectric applications 1 (power gen- 
eration or cooling). The recent focus has been on inves- 
tigating strongly correlated materials, which may prove 
to have better performance than the current generation 
of semiconductor-based devices. The two main areas of 
application for thermoelectrics are in power generation 
from the Peltier effect, 2 where heat energy is converted 
into electricity; and in thermoelectric cooling, where an 
electrical current is driven through a device to force heat 
to move from the cold to hot end. Power generation ap- 
plications typically operate at temperatures higher than 
1000 K, with the heat source being a radioactive material 
(for applications in the space industry) or a combustion 
source. Thermoelectric coolers usually operate around 
room temperature, and the semiconductor-based devices 
do not function below about 200 K. Currently, thermo- 
electric devices fit niche markets, where reliability, size, 
or weight are more important than efficiency. The cool- 
ers usually operate with relatively low heat loads because 
of their poor efficiency. 

The efficiency of a thermoelectric device is a function 
of the dimensionless product of a materials parameter 
denoted Z with the average temperature T (between the 
hot and cold heat sources of the device) and is called ZT. 
It satisfies 



ZT = 



Ta dc S 2 

K e + Hi 



(1) 



and the term ad c S 2 in the numerator is often called the 
power factor. Here we have Udc the dc electrical conduc- 
tivity, S is the Seebeck coefficient, 3 T is the temperature, 
K e is the electronic contribution to the thermal conduc- 
tivity, and Ki is the lattice contribution to the thermal 
conductivity (we are assuming the electron-phonon inter- 
action is small enough that these two effects can be de- 
coupled). It is commonly stated that ZT > 1 is needed 
for operation of thermoelectric devices, but this is not 



necessarily true. 4 For example, if we consider a thermo- 
electric cooler operating at 300 K and with 50 K of cool- 
ing, then one can operate such a device for ZT > 0.7, 
but one would need ZT ss 4 to achieve the coefficient of 
performance of a conventional compressor-based refrig- 
erator (which lies in the 1.2 — 1.4 range). Nevertheless, 
most commercial thermoelectric devices have ZT near 
1 because no materials have been discovered with much 
larger values. Of course there is significant interest in 
increasing ZT by a factor of 4 at room temperature (to 
be competitive with conventional coolant-based technol- 
ogy), or to above 1 at low temperature to allow for new 
applications such as solid-state coolers for superconduct- 
ing electronics or infrared detectors. 

Although there are no fundamental thermodynamic 
limitations 5 to the size of ZT, it has proved to be quite 
difficult to find materials with ZT > 1 over a wide tem- 
perature range, and to find much larger values of ZT 
(say ZT > 3). Recently, Rontani and Sham 6 proposed 
that heterostructures of correlated semiconductors and 
metals could have dramatically large values of ZT at low 
temperature. Their idea was that if one tuned the large 
electronic density of states (DOS) of the /-electrons to lie 
close to the fermi level, then one could produce huge val- 
ues of ZT (earlier work proposed similar ideas as well 7 ). 
Mahan and Sofo 8 also argued in the same vein for opti- 
mization in bulk materials in 1996. But so far, no one has 
been able to demonstrate that such large values of ZT 
are possible in a true many-body system (and are not the 
artifact of some approximations employed in the analy- 
sis). We examine this problem in detail for the Falicov- 
Kimball model here. By working in the limit where the 
ion concentration is fixed and nonzero at T — 0, we have 
adjusted the renormalized energy level of the ion to lie at 
the electronic chemical potential, which has the potential 
for producing large thermoelectric responses. 

Mahan and Sofo 8 also proved that the figure-of-merit 
always satisfies an inequality 



ZT < 



Ta dc S 2 



Kl 



(2) 



1 



regardless of the strength of the many-body interactions. 
This has important implications for theorists, because in 
purely electronic models, such as the one we investigate 
here, m = 0, so there is no a priori limitation on the 
magnitude of ZT. But it also presents some problems 
for low-T calculations, since the electronic contribution 
to the thermal conductivity is usually much smaller than 
the lattice contribution at low temperature (especially 
for insulators), and hence purely electronic estimates of 
ZT can be greatly enhanced when the lattice effects are 
ignored. This becomes less of an issue at high temper- 
ature, where the electronic contribution to the thermal 
conductivity can dominate. 

Another interesting feature that plagues the low-T 
regime is the fact that in most systems S — ► as T — > 0. 
Since the ratio of the conductivities often satisfies the 
Wiedemann-Franz law 



H 



Ode 



k B 

e 



CT 



(3) 



with C the Lorenz number (equal to 7r 2 /3 in a fermi liquid 
and 3 in an intrinsic semiconductor) , we have ZT — > if 
S —* at low temperature. Similarly, if S(T) suffers a 
sign change at any T, then ZT will be quite low in the 
vicinity of the sign change. 

The Falicov-Kimball model 9 appears to be able to de- 
scribe an increasing number of materials and systems. 
One example, that fits within the binary-alloy picture, is 
tantalum dcfficicnt tantalum nitride 10 Ta 2 N. This mate- 
rial is metallic when x = 1 but becomes a fairly large-gap 
insulator (about 1.5 ev) when x = 0.6. If we let the A ion 
denote a unit cell with a Ta atom, and a B ion denote 
a unit cell with no Ta atom, then U is the difference in 
site energies for the two configurations. The total con- 
duction electron concentration also depends on the Ta 
vacancies, as each vacancy can bind 5 electrons. It is 
easy to model the metal-insulator transition at x = 0.6 
by properly varying the electron concentration with the 
concentration of Ta vacancies. 

In Section II we develop the formalism for deriving 
the dc conductivity, the thermopower and the thermal 
conductivity. We use this to determine both the Lorenz 
number and the figure-of- merit. In Section III we provide 
numerical results for the thermal transport illustrating 
regimes where ZT can become large, and describing the 
physical principles that drive such enhancements. In ad- 
dition, we describe in detail the situation behind a large 
figure-of-merit at low-T and whether achieving such a 
goal is feasible. Conclusions are presented in Section IV. 



II. FORMALISM FOR THE THERMAL 
TRANSPORT 

The Hamiltonian we study is the spinless Falicov- 
Kimball model 9 with a canonical-binary-alloy picture 



(hj) 2 



(4) 



where c\ (cj) is the electron creation (annihilation) oper- 
ator for a spinless electron at site i (spin can be included 
trivially if desired by doubling all of the defined be- 
low), Wi is a variable that equals zero or one and cor- 
responds to the presence of an A ion (wi = 1) or the 
presence of a B ion (wi = 0) at site-i, and U is the in- 
teraction strength (difference in the site energy between 
the A and B ions). The hopping integral is scaled with 
the spatial dimension d so as to have a finite result in the 
limit 11 d — > oo; we measure all energies in units oft* = 1. 
We work on a hypercubic lattice where the noninteracting 
density of states is a Gaussian p(e) = cxp(— e 2 )/ y/irt*fl uc 
(with Q uc the volume of the unit cell). A chemical po- 
tential /i is employed to adjust the conduction electron 
filling p e . 

The Falicov-Kimball model can be solved exactly by 
employing dynamical mean field theory 12,13 . A review 
that describes how to solve for the Green's function 
using the equation of motion technique has recently 
appeared. 14 Because the self energy E(z) has no momen- 
tum dependence, the local Green's function satisfies 



G(z) = J dep(e) 



1 



z + fi - E(z) 



(5) 



with z a complex variable in the complex plane. The self 
energy, local Green's function, and effective medium Go 
are related to each other by 



G - 1 (z)-G- 1 (z) = S(z), 
and the Green's function also satisfies 
G(z) = (1 — w\)Gq(z) + wi— j- 1 



G^(z)-U 



(6) 



(7) 



Here w\ is the average concentration of the A ions (which 
is an input parameter). The algorithm for determining 
the Green's function begins with the self energy set equal 
to zero. Then Eq. (5) is used to find the local Green's 
function. The effective medium is found from Eq. (6). 
The new local Green's function is then found from Eq. (7) 
and the new self energy from Eq. (6). This algorithm is 
repeated until it converges. 

When these equations are solved we find a number 
of interesting results for the single-particle properties. 
First, both the interacting DOS [pmi(w) = — Im{G(o; + 
z<5)} / 7r] and the self energy on the real axis are indepen- 
dent of temperature 15 when w\ and /x are fixed (all of the 
temperature dependence for fixed p e arises from the tem- 
perature dependence of /i, which shifts the zero- frequency 
point of the DOS). Second, we find that the self energy 
does not display fermi-liquid properties unless [7 = 0, 
Wi = or wi = 1 (which are all noninteracting cases). 
In particular, we do find (for small enough U) that the 
imaginary part of the self energy is quadratic in lo, but 



2 



the curvature has the wrong sign, and the zero frequency 
value of the imaginary part of the self energy does not 
go to zero as T — > (in fact, it remains fixed for all T). 
Third, the real part of the self energy is linear (for small 
enough U), but the slope has the opposite sign of what 
is seen in a fermi liquid. Finally, we see that if U is large 
enough, the DOS develops a gap, and if the electronic 
chemical potential lies in the gap, then the self energy has 
quite anomalous behavior (including developing a pole). 

Transport properties are calculated within a Kubo- 
Greenwood formalism 16 . This relates the transport co- 
efficients to correlation functions of the corresponding 
transport current operators (these are equal to the bare 
bubbles because there are no vertex corrections in the 
large-dimensional limit 17 ). We need two current opera- 
tors here — the particle current 18 



J 



E 



(8) 



(where the velocity operator is v g = V q e(q), the band- 
structure is e(q), and the Fourier transform of the cre- 
ation operator is cj = . exp[i<? • Rj]cJ/N) and the heat 
current 19 ' 18 



(9) 



[where W(q) = Ylj exp(— • TLj)wj/N]. 

The dc conductivity, thcrmopowcr and electronic ther- 
mal conductivity can all be determined from relevant 
correlation functions of the current operators. We de- 
fine three transport coefficients Ln, L 12 = L 2 i, and L 2 2- 
Then 



2 T 



S = 



k R L 



12 



lelTLn' 



and 



h 2 



L 



L\ 2 L 



21 



22 



41 



(10) 



(11) 



(12) 



One finds that the electric and thermal conductivities are 
always positive, but the thermopower can have either 
sign — a positive thermopower corresponds to electron- 
like transport and a negative thermopower to hole-like 
transport (we use the sign convention of Ashcroft and 
Mermin 20 ). The transport coefficients are found from 
the analytic continuation of the relevant "polarization 
operators" at zero frequency 



L\\ = lim Re— Ln(i>), 

v—>0 V 



Lu(ivi) 



f 

Jo 



dre^(T T ji( T )j fi {0)), 



(13) 



where v\ = 2ttTI is the Bosonic Matsubara frequency 
{fi = 1/T), the r-dependence of the operator is with re- 
spect to the full Hamiltonian in Eq. (4), and we must 
analytically continue Ln(ivi) to the real axis L\\(v) be- 
fore taking the limit v — * 0. Similar definitions hold for 
the other transport coefficients: 



L12 = lim Re-L^M, 

v— >o v 



L\ 2 (in) 



f 

Jo 



dre^(T T ji(T)j Qf3 (0)), 



(14) 



and 



L 22 = lim Re-L 2 2(^), 

v— >0 V 

U 2 (iv x ) = f dre^{T T ^ Qa {r)j Q p{Q)). (15) 
Jo 

In all of these equations, the subscripts a and (3 denote 
the respective spatial index of the current vectors. Note 
that our definitions of the Lij coefficients, while standard, 
have one less factor of T than seen in some other papers. 

The analytic continuation is straightforward, and pro- 
duces the Mott form 21 for the transport coefficients 
(which is usually called the Jonson-Mahan theorem, 19 
and was explicitly evaluated for the Falicov-Kimball 
model 22 ): 



- _ao r 
H3 e 2 L 



dto 



dm 

dui 



(16) 



with the relaxation time r(w) defined by 
ImG(w) 



t(u) 



+ 2-2Rc{[w + ^-S(w)]G(a;)} (17) 



f(uj) = 1/[1 + exp(/3w)] and <r = e 2 ir 2 /hda d ~ 2 on a hy- 
percubic lattice in d-dimensions. Note that even though 
we represented the above form by an effective relaxation 
time, the above expression for the transport coefficients 
is exact. 

Once the transport coefficients are known, we can de- 
termine the electrical and thermal conductivities and the 
thermopower. In addition, we find the figure-of-merit 
satisfies 



ZT = 



r 2 

^12 



r 2 



LuL 22 

when we neglect ki and the Lorenz number becomes 



(18) 



' e) -a dc T-{e) L\ X T 2 ' [W) 



Note that because —df(u>)/du> is an even function of 
u), only the even part of r(w), t £ (uj) = \t{uj) +t(— w)]/2, 



3 



contributes to the Ln and L22 coefficients, and only the 
odd part t (u>) = [t(u>) — t(—u>)]/2 contributes to Li 2 . 
Furthermore, one expects on general physical grounds 
that t(cj) > 0. Indeed, this can be easily shown to be 
true by rearranging Eq. (17) into 



2[ImS(w) 



{[lu + p - ReE(w) - e] s 



[ImS(w)] 2 } 2 
(20) 



which is manifestly nonnegative. In addition, we see 
that t(lu) vanishes whenever Im£(w) = 0. In the case 
of particle-hole symmetry, when p e = w\ — 0.5, it is easy 
to show that t {lo) — 0, and both the thermopower and 
ZT must vanish. 

We can employ the nonnegativity of r(w) to show the 
Mahan-Sofo 8 bound for ZT. The starting point is the 
fact that 



3/ <M- 



dui 



t(w)(u> + a) 2 > 



(21) 



which holds since all terms in the integrand are positive. 
Choosing a = Li 2 /Ln then yields LuL 2 2 — L\ 2 > 
(equality only occurs if the integrand is a delta function, 
which has no variance). Now the full expression for ZT 
can be written as 



ZT = 



T 2 
^12 



/L11L2 



1 - L? 2 /(LnL22) + /s,T/(fe|L 22 ) 1 



(22) 



by dividing the numerator and denominator in Eq. (1) 
by k%L22/T. Here we used the notation £ = L 2 2 /LuL 2 2 
and A = k{I '/ '{k 2 B L 22 ). Using the fact that £ < 1 then 
tells us that £A< A<l-£ + A since l-£ > 0. Dividing 
both sides of the inequality by A(l — £ + A) then yields 



ZT - 



1 

— < — 

A - A 



n e + Ta dc S 2 



(23) 



which is the Mahan-Sofo result. The inequality in 
Eq. (23) holds for any system that satisfies the Jonson- 
Mahan theorem, regardless of the many-body interac- 
tions present, as long as there is no electron-phonon cou- 
pling, which precludes the separation of the thermal con- 
ductivity into electronic and lattice pieces. 

The integrals for the transport coefficients all have a 
derivative of the fermi function in them. This derivative 
becomes strongly peaked around u> — with a width 
on the order of T as T — > 0. In metals, we typically 
find that the relaxation time can be written in the form 
t{u>) = tq + wt 1 + O{oj 2 ) for small to. If we include only 
the first two terms in the expansion for the relaxation 
time, we find that 

„ k B ttVT 



k% n 2 cr T T 



(24) 



to lowest-order in T. We can read off that the Lorenz 
number C is equal to 7r 2 /3 and ZT — ► tt 2 t' 2 T 2 /3 (when 
we neglect «;). Hence, ZT will be small in metals 
at low temperature because Wiedemann-Franz law fixes 
the ratio of conductivities, and the thermopower van- 
ishes as T — > 0. If the Wiedemann-Franz rule holds 
at higher temperature as well, then one needs to find 
a thermopower larger than nks /V3\e\ which is equal to 
156 /tiV/K in order to have ZT > 1 in a metal. There 
are no known metals that have thermopowers larger than 
125 /xV/K and most metals are one to two orders of mag- 
nitude smaller, 5 so the only way to find metals that are 
useful for thermoelectric power generation or cooling is 
if the Wiedemann-Franz law does not hold. 

The situation in an insulator is quite different. If we 
assume that t{lo) = Topi n t(^) with to a constant, and 
if we choose a generic interacting density of states that 
increases like a power law, 

p int {u) = e{io - E g /2)C(Lu - E g /2) a 

+ 6(-Lu-E g /2)C'(-oj-E g /2) a ' (25) 

[with 9{x) the unit step function and E g the insulating 
gap], then the behavior of the transport coefficients dif- 
fers from Eq. (24) . The starting point is to note that the 
number of holes excited in the lower band is equal to the 
number of electrons in the upper band, and at low tem- 
perature the fermi factors can be replaced by Boltzmann 
factors, resulting in 



diop int (cj)e^- A ^ 



(Lop int (w)e-^- A ^ 



(26) 



with the chemical potential written as p, = po + Ap(T). 
Solving for Ap(T) yields 



Ap(T) 



/ dojp int (Lu)exp((3u) 
.J —00 

dujp in t{io) exp(-/3w) 



(27) 



with po the chemical potential at T = 0, which lies in 
the middle of the insulating gap (and is chosen to be the 
origin here). Using the form for the DOS in Eq. (25), 
then yields 



M?1 = — TlnT+-ln 



C'T(a' + 1) 
CT(a + 1) 



(28) 



This form for Ap(T) plays an important role in the 
thermal transport of insulators. Indeed, an examination 
of the in and Li 2 coefficients in an insulator, with T <SC 
E g /2, leads immediately to 



£12 = Li 



Ap + p 



dAp 



(29) 



4 



which follows by taking the derivative of Eq. (26) with 
respect to (3 and using the definitions for L\\ and L\i 
in the low-temperature limit. If A/i(T) = 0, as occurs 
at half filling, then we immediately see that S = (as 
we also know must hold due to particle- hole symmetry). 
Furthermore, if Ap(T) depends linearly on T, the ther- 
mopower vanishes linearly in T as well. S can approach a 
constant value if the low-temperature slope diverges [like 
in the TlnT dependence in Eq. (28)]. Divergent behav- 
ior can occur if A/i(T) depends on T with a power-law 
less than 1 (and greater than 0), but it seems unlikely 
that one could find a circumstance where S will diverge 
like 1/T at low temperature (note that this argument 
does not hold for disordered insulators, that have local- 
ized states, so that the DOS does not vanish at the fermi 
level 23 ). Hence the qualitative temperature dependence 
of the thermopower in a correlated insulator can be de- 
termined by the temperature dependence of the chemical 
potential (at low T). Note that the above arguments ob- 
viously hold for intrinsic semiconductors as well. 

If we use the generic insulator DOS in Eq. (25) and 
the approximation that t(uj) — T p int (uj), then we find 
the following behavior for the transport in an insulator 



_<7qto 0Ea/2 C'TTja' + l) 
dc 4 V CT<*T{a+l) 

S = ^r(a-a') + 0(T) 



K _ Wo p -0E g /2 C'T°'T(a> + l) 
4 V CT<*T{a + l) 



k 2 



E 2 

-1 - (a 1 - a)T 
2T y ' 



(30) 



The Lorenz number then diverges like C = E 2 /2T 2 at 
low temperature. 

In a correlated insulator, we must have a — a', so the 
thermopower vanishes linearly as T — > 0. This occurs be- 
cause the interacting density of states develops a gap due 
to a pole in the self energy. Hence the band edge should 
have the same power law as that of the noninteracting 
system, which is determined by the van Hove singularity, 
and must be the same for the upper and lower band edge; 
the only place where this analysis might break down is 
at the critical coupling strength, where the gap initially 
opens, but we expect that the power-law dependence of 
the band edge may differ from the generic case, but it 
still should be equal for the upper and lower branches. 

Unfortunately the hypercubic lattice in large dimen- 
sions does not behave like the generic insulator in the 
strong-coupling limit. This happens because the DOS 
never vanishes for a finite range of lo, but rather ap- 
proaches zero exponentially for all nonzero lo (it is sup- 
pressed to zero precisely at lo = and T = 0). Hence we 
need to analyze the situation for the Gaussian DOS in 
more detail. 

Inside the "gap region" of the DOS, the real part of 



the self energy is large, because the self energy develops 
a pole at lo = in a correlated insulator. Similarly, the 
imaginary part of the self energy is very small. Hence, we 
can approximately determine the Hilbert transformation 
in Eq. (5) as 



G( W ) 



inp[oj + /i- Re£(u;)]. (31) 



to + /j, — ReE(w) 

Then we find from Eq. (6) that 

Gq 1 {lo) ~ lo + ji + i{ImE(oj) 

+ 7r[w + ^-RcE(w)] 2 / 9[^ + A*-RcS(w)]}. (32) 

Next, by employing Eqs. (7) and (6), we finally determine 
that 



ReS(w) 



W\U{bJ + /l) 

lo + p — (1 — wi)U 



(33) 



and 

ImS(w) 



(uj + p) 2 (uj + fi-U) 2 
' [lu + /i - (1 - Wl )U} 2 P 
wi(l - wi)U 2 



(u) + /i)(w + /i - U) 

LO + n — (1 — Wi)U 



[lo + /j, - (1 - Wl )Uf + wi(l - w^U 2 ' 



Now the local self energy of the correlated insulator has 
a pole at lo = 0, hence we learn that = (1 — wi)U 
for the correlated insulator, and the interacting DOS is 
equal to zero at lo = (and T = 0). To lowest order, we 
find that the scattering time then satisfies 



-{lo) 



[lo + n - (i - Wl )uy 



{lo + h) 2 (lj + ii - C/) 2 wi(l - Wl )U 2 



(35) 



which is a quartic dependence on frequency in the "gap 
region." This implies that for the hypercubic lattice, the 
scattering time is much bigger than what one would guess 
from the approximation t{lo) = Top int {uj) (which would 
be exponentially small). 

We find that in our calculations, we usually have to 
input the imaginary part of the Green's function and of 
the self energy by hand in the "gap region" , because the 
numerical precision of the computer does not allow them 
to be determined accurately with the iterative algorithm 
(the difficult step is constructing the self energy from the 
difference of the inverse of the local Green's function and 
the effective medium, where numerical precision is lost). 
For accurate calculations, we need to go to higher order 
in our expansion than what was illustrated above, and 
since the real part of the self energy is determined reliably 
in the iterative algorithm, we construct the imaginary 
parts of the Green's function and the self energy from the 
actual converged value of the real part of the self energy, 
rather than employing the asymptotic forms discussed 
above. Even with all of these tricks, it still is a challenge 
to determine a smooth t{lo) for the correlated insulator. 

There is an advantage to studying the Falicov-Kimball 
model over the more popular Hubbard model, 24 because 



(34) 



5 



the Falicov-Kimball model can be tuned to have a metal- 
insulator transition for any value of w\ , simply by choos- 
ing U large enough and fixing p e = l—w\ [p, = (1— wi)U\. 
The Hubbard model has a metal-insulator transition only 
at half filling, where the thermopower vanishes due to 
particle-hole symmetry (this can be broken by introduc- 
ing multiband Hubbard models 25 ). Hence we can study 
effects of the correlated metal-insulator transition in the 
Falicov-Kimball model, that are inaccessible in the single- 
band Hubbard model. Since real materials typically have 
complicated bandstructures, one does not expect them to 
be particle-hole symmetric except in very special circum- 
stances. Once again, the Falicov-Kimball model can be 
viewed as a more generic metal-insulator transition for 
making contact with real materials. The only disadvan- 
tage is that the Falicov-Kimball model is a non-fcrmi- 
liquid except in "noninteracting" limits (where it is a 
fermi gas). 



III. NUMERICAL RESULTS 

We present results at p e = 1 — w\ for three different 
values of U: (i) U = 1 which is a strongly correlated 
metal, that has a dip or kink in the interacting DOS at 
the fermi level; (ii) U = 1.5 which undergoes a metal- 
insulator transition as a function of w\\ and (iii) U = 2 
which is a correlated insulator with a sizable "gap re- 
gion" . 
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Frequency gj [t*] 

FIG. 1. (a) Effective scattering time and (b) interacting 
DOS for the Falicov-Kimball model at U = 1, p e = 1 — wi, 
and wi = 0.5 (solid), 0.4 (dashed), 0.3 (chain-dotted), 0.2 
(dotted), and 0.1 (chain-triple-dotted). Note how both t(uj) 
and the DOS develop dips at the fermi level; both functions 
are symmetric at half filling uii = 0.5. These plots have fixed 
the origin at the T — > limit of the chemical potential. 



The effective scattering time and the interacting DOS 
are plotted in Fig. 1. We choose U = 1, p e = 1 — w\, and 
vary w\ (0.5, 0.4, 0.3, 0.2, and 0.1). Note how the strong 
correlations create a dip in both the interacting DOS and 
the relaxation time near the fermi level. As w\ is made 
smaller, the scattering time becomes more asymmetric in 
frequency, which should yield larger thermopowers. Note 
that there is residual scattering at the fermi energy as 
cij — > 0, which should produce a finite value for the DC 
conductivity at T = (recall the Falicov-Kimball model 
is not a fermi liquid in this regime). 
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FIG. 2. (a) DC conductivity (in units of cro), (b) ther- 
mopower (in units of fes/|e|), and (c) electronic contribution 
to the thermal conductivity (in units of k%a /e 2 ) for the Fal- 
icov-Kimball model at U = 1, p e = 1 — wi, and wi = 0.5 
(solid), 0.4 (dashed), 0.3 (chain-dotted), 0.2 (dotted), and 0.1 
(chain-triple-dotted) . 

We plot the DC conductivity, thermopower, and elec- 
tronic contribution to the thermal conductivity in Fig. 2. 
The DC conductivity is fairly flat as a function of T over 
a wide temperature range, and is enhanced, as w\ is 
made smaller (and the particle-hole asymmetry is en- 
hanced). The thermopower is essentially flat at high 
temperature, increases as w\ decreases, and then has 
a crossover to linear behavior below T 0.25i*. The 
slope of the low-temperature linear behavior increases 
as particle-hole asymmetry increases. The thermal con- 
ductivity is also essentially flat at high temperature, but 
then has a characteristic linear behavior that sets in be- 
low T w 0.6i*. This behavior is typical of a metal that 
has significant scattering. Note that at low temperature, 
we might expect the Wiedemann-Franz law to hold, while 
at high temperature, the Lorenz number must decrease 
with T (since the ratio of the conductivities approaches 
a constant). The high temperature decrease of C with a 
constant (or increasing) thermopower can lead to a large 
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FIG. 3. (a) Lorenz number and (b) thermoelectric figure of 
merit for the Falicov-Kimball model at U = 1, p e = 1 — wi, 
and wi = 0.5 (solid), 0.4 (dashed), 0.3 (chain-dotted), 0.2 
(dotted), and 0.1 (chain-triple-dotted). 

The Lorenz number and thermoelectric figure of merit 
are plotted in Fig. 3. Note how the Lorenz number ap- 
proaches 7r 2 /3 as T — > 0. This follows whenever t(uj) ^ 
at u = 0, or whenever the interacting DOS at the chem- 
ical potential is nonzero. The decrease in C at high T 
occurs because both the DC conductivity and the elec- 
tronic contribution to the thermal conductivity become 
flat at high T. What is interesting is the moderate tem- 
perature peak in C with the linear decrease toward the 
fcrmi-liquid value as T — > 0. This occurs because the 
system has strong scattering, which produces deviations 
from the Wiedemann- Franz law at moderate T. What is 
unfortunate is that an enhancement of C leads to a re- 
duction in ZT. This can be seen clearly in panel (b). The 
high-temperature figure of merit is large, and increases 
in magnitude as w\ decreases due to the enhancement of 
the thermopower (since C is independent of w\ at high 
T). But the figure of merit rapidly decreases at low-T 
and indicates that strongly scattering metals of this type 
can not be used for thermoelectric applications at low 
temperature (but they would work fine at high T if kl 
was small enough, which is usually true). 
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FIG. 4. (a) Effective scattering time and (b) interacting 
DOS for the Falicov-Kimball model at U = 1.5, p e = 1 — wi, 
and wi = 0.5 (solid), 0.4 (dashed), 0.3 (chain-dotted), 0.2 
(dotted), and 0.1 (chain-triple-dotted). Now t(ui) has a quar- 
tic dependence on frequency, and the DOS becomes exponen- 
tially small in the "gap region". The gap fills in for wi < 0.25 
(which is difficult to see in the figure), so we expect the trans- 
port behavior to be different for small wi . 

Next we show the evolution of the relaxation time and 
the interacting DOS as we move into the correlated in- 
sulator regime U = 1.5. Here the self energy develops a 
pole at oj = for w\ > 0.25, which produces a narrow 
region of exponentially small DOS inside a "gap region" , 
but the pole disappears for small enough w\, and the sys- 
tem has strong scattering, but is still metallic as T — > 0. 
This is difficult to see in the above figure because we are 
not plotting the gap region on a logarithmic scale. We 
expect that the transport properties may differ depend- 
ing on the value of w\ for U — 1.5, but it should not be 
too dramatic because the metallic phase has very strong 
scattering. 
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FIG. 5. (a) DC conductivity, (b) thermopower, and (c) 
electronic contribution to the thermal conductivity for the 
Falicov-Kimball model at U = 1.5, p e = 1 — wi, and wi = 0.5 
(solid), 0.4 (dashed), 0.3 (chain-dotted), 0.2 (dotted), and 0.1 
(chain-triple-dotted). Note the very low-temperature peak in 
the thermopower for t«i = 0.4 and wi = 0.3. 

The transport coefficients are plotted in Fig. 5. These 
look similar to those in Fig. 2 except now the conductiv- 
ity vanishes as T — > for the insulating phases and the 
thermal conductivity approaches zero faster than linearly 
for the insulators as well. In addition, when w\ is close 
to half filling, we find a low-T peak in the thermopower. 
Our calculations indicate that the thermopower goes to 
zero as T — > 0, but we cannot rule out the possibility of a 
small nonzero value at T = (the chemical potential, as 
a function of T, develops a large negative slope as T — > 0, 
but we cannot tell if the slope diverges). 
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FIG. 6. (a) Lorenz number and (b) thermoelectric figure of 
merit for the Falicov-Kimball model at U = 1.5, p e = 1 — wi, 
and wi = 0.5 (solid), 0.4 (dashed), 0.3 (chain-dotted), 0.2 
(dotted), and 0.1 (chain-triple-dotted). 

The plot of the Lorenz number and of the thermoelec- 
tric figure of merit in Fig. 6 show a number of interesting 
features. First off, note the very large increase in C as 
T — ► for the half-filled case versus all other fillings. This 
behavior happens because the half-filled case is qualita- 
tively different from all other cases at low temperature. 
At half filling, the chemical potential is always at U/2 
and has no temperature dependence. Hence the system 
always has a pseudogap-likc behavior of the DOS, with 
it vanishing only at the fermi level. For all other cases, 
the chemical potential depends on T. Hence, at finite- 
ly the DOS at the chemical potential is nonzero, and 
only vanishes exactly at T = 0. This makes the low- 
temperature behavior away from half filling appear like 
that of a metal, resulting in the Wiedemann-Franz pre- 
diction for C. Hence there is a low-T downturn to C off 
of half filling. Second, at high temperature, C has dop- 
ing dependence now, which is a result of the fact that 
for W\ < 0.25 the system is metallic, while for wi > 0.25 
it is insulating. Third, the values of the thermoelectric 
figure of merit are enhanced at moderate temperature 
for low tui, and there is a low-temperature peak in ZT 
for w\ = 0.3 and 0.4, which is hard to see in the figure 
(it is blown up in Fig. 10). Finally, the high-temperature 
behavior of ZT is actually reduced over that of the corre- 
lated metal because C is somewhat larger and S is some- 
what smaller here. ZT is also a flatter function of T 
here. 
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FIG. 7. (a) Effective scattering time and (b) interacting 
DOS for the Falicov-Kimball model at U = 2, p e = 1 — wi, 
and wi = 0.5 (solid), 0.4 (dashed), 0.3 (chain-dotted), 0.2 
(dotted), and 0.1 (chain-triple-dotted). 

We plot the relaxation time and interacting DOS for 
the case U = 2 in Fig. 7. Now all dopings have well- 
defined "gap regions" , but the size of the "gap" is still 
relatively small. Note that the effective relaxation time 
is quartic for all fillings here. 
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FIG. 8. (a) DC conductivity, (b) thermopower, and (c) 
electronic contribution to the thermal conductivity for the 
Falicov-Kimball model at U = 2, p e = 1 — wi, and wi = 0.5 
(solid), 0.4 (dashed), 0.3 (chain-dotted), 0.2 (dotted), and 0.1 
(chain-triple-dotted) . 



The transport coefficients are plotted in Fig. 8 for 
[7 = 2. Both the electrical and thermal conductivities be- 
have as expected, with exponentially small values at low 
T (but the "gap" decreases as w\ decreases, so the ex- 
ponent is doping dependent). The thermopower still has 
an interesting low-temperature peak. The peak height 
moves out to larger T and broadens as U is increased 
and as w\ is decreased. We expect this will have a sig- 
nificant impact on ZT . 
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FIG. 9. (a) Lorenz number and (b) thermoelectric figure of 
merit for the Falicov-Kimball model at U = 2, p e = 1 — wi, 
and wi = 0.5 (solid), 0.4 (dashed), 0.3 (chain-dotted), 0.2 
(dotted), and 0.1 (chain-triple-dotted). 

The Lorenz number becomes huge (but does not ap- 
pear to diverge) in the correlated insulator at half filling, 
but continues to have the low-T peak and return to 7r 2 /3 
as T — > away from half filling. The thermoelectric fig- 
ure of merit has a sharp peak at low temperature, whose 
width and peak location increase with an increase of U 
and a decrease of w\. Note that the low-T peak is as- 
sociated with the sharp drop in C at low temperature, 
which seems to occur only when one has pseudogap be- 
havior and a chemical potential that moves away from 
the point where the DOS vanishes at finite T. Never- 
theless, the peak is rather striking and does show that 
ZT > 1 is possible in an all electronic system. Since the 
electronic contribution to the thermal conductivity expo- 
nentially decreases at low temperature, while the lattice 
contribution decreases as a cubic power law (when one 
is well below the Debye energy), we expect that at very 
low temperatures the thermal conductivity will be domi- 
nated by a lattice contribution, and if that contribution is 
significantly larger than the electronic contribution, then 
the low-T peak in ZT will go away. 
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FIG. 10. Blow-up of the low-temperature region for the 
thermoelectric figure of merit with (a) U = 1.5 and (b) {7 = 2. 
The curves correspond to wi = 0.5 (solid), 0.4 (dashed), 
0.3 (chain-dotted), 0.2 (dotted), and 0.1 (chain-triple-dotted). 
Note how the peak can be pushed to very low T and to rea- 
sonably high values when U tuned to lie close to the critical U 
of the metal-insulator transition. The peak generically moves 
to higher T, increases in magnitude, and broadens as U is 
increased. 

We blow up the ZT plots for the low temperature re- 
gion in Fig. 10. Note how one can find a sharp peak with 
an enhanced ZT at low temperature when one is close 
to the metal-insulator transition. As U is increased, the 
magnitude of the peak increases and broadens, but it is 
pushed to higher values of temperature. We expect that 
the lattice contribution to the thermal conductivity may 
significantly reduce this peak, or may even destroy it, 
but in both cases, the thermal conductivity at the low-T 
peak is about 10~ 4 the value at high T and perhaps the 
lattice thermal conductivity could be low enough that it 
would not significantly interfere with the peak. 
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FIG. 11. (a) DC conductivity, (b) thermopower, and (c) 
electronic contribution to the thermal conductivity for the 
Falicov-Kimball model at U = 1, wi = 0.2, and p e + wi = 0.8 
(solid), 0.9 (dashed), 1.0 (chain-dotted), 1.05 (dotted), and 
1.1 (chain-triple-dotted). 

We conclude our presentation of results by showing 
what happens when the system is forced to be metallic, 
by pushing the electronic chemical potential away from 
the "gap region" . We start with the strongly correlated 
metal at U = 1 and w\ = 0.2, but we vary the total 
concentration to be 0.8 (solid), 0.9 (dashed), 1.0 (chain- 
dotted), 1.05 (dotted), and 1.1 (chain-triple-dotted). The 
relaxation time and the interacting density of states can 
be read off of the w\ = 0.2 curves in Fig. 1, with the only 
change an overall shift of the origin for the different elec- 
tronic fillings. The transport coefficients are plotted in 
Fig. 11. As expected, the system is least conductive when 
the filling is equal to 1 (corresponding to the chemical po- 
tential at the "kink" in the interacting DOS). The ther- 
mopower is enhanced at low temperature as we approach 
a total filling of one from below and suppressed at high 
temperature. The converse occurs for fillings larger than 
1 — the low-temperature thermopower is reduced and the 
high-temperature is enhanced. This suggests that the 
thermoelectric properties at high temperature may be 
enhanced by increasing the electron filling. 
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FIG. 12. (a) Lorenz number and (b) thermoelectric figure 
of merit for the Falicov-Kimball model at U = 1, wi = 0.2, 
and p e + wi = 0.8 (solid), 0.9 (dashed), 1.0 (chain-dotted), 
1.05 (dotted), and 1.1 (chain-triple-dotted). 

The Lorenz number and thermoelectric figure of merit 
are plotted in Fig. 12. Note how the peak in the Lorenz 
number only occurs for total fillings larger than about 
0.85, and how it gets larger and moves to higher tem- 
perature as the electronic filling increases. This has an 
obvious effect on ZT: we find ZT is enhanced at low 
temperature for fillings less than 1, and is enhanced at 
high temperature for fillings greater than 1. These effects 
can be significant. 
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FIG. 13. (a) DC conductivity, (b) thermopower, and (c) 
electronic contribution to the thermal conductivity for the 
Falicov-Kimball model at U = 3, wi — 0.2, and p e + wi = 0.8 
(solid), 0.9 (dashed), 1.0 (chain-dotted), 1.05 (dotted), and 
1.1 (chain-triple-dotted). 

Finally, we examine a strongly correlated insulator 
with U = 3 and W\ = 0.2. We don't plot the relax- 
ation time and DOS here, because it is similar to what 
is seen for U — 2, but with a somewhat larger "gap re- 
gion". The transport coefficients are plotted in Fig. 13. 
Here we see a marked difference between the metallic 
cases, with total filling not equal to 1 and the insulating 
case where it equals 1. In particular, the DC conductiv- 
ity rises as T decreases, as expected for a metal, but it 
is exponentially suppressed in the insulator. The ther- 
mal conductivity is similar, with well-developed linear 
regimes for the metallic systems (that are better defined 
for fillings less than 1), and "gapped" behavior for the in- 
sulator. The thermopower has quite interesting behavior. 
For fillings less than 1, we see a dramatic enhancement 
in the low-temperature thermopower as we approach the 
insulator. For fillings above 1, we see the thermopower 
has a sign change. This sign change is easy to under- 
stand. When the total filling lies in the range between 
1 and 1.1, the electronic chemical potential lies in the 
lower half of the upper Hubbard band at low tempera- 
ture. Hence at low-T, the thermopower should appear to 
be electron- like (positive). But as T is increased, the in- 
fluence of the gap becomes smaller, and the system looks 
overall hole-like because the chemical potential is in the 
top part of the overall "bandstructure" . The crossover 
temperature should be on the order of the size of the 
"gap region". In addition to this sign change, we see a 
large enhancement of the high-temperature thermopower 
as the filling increases. This is also expected because the 
chemical potential is moving higher and higher in the 
band. 
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FIG. 14. (a) Lorenz number and (b) thermoelectric figure 
of merit for the Falicov-Kimball model at U — 3, wi = 0.2, 
and p e + wi = 0.8 (solid), 0.9 (dashed), 1.0 (chain-dotted), 
1.05 (dotted), and 1.1 (chain-triple-dotted). 

We end with the Lorenz number and thermoelectric 
figure of merit for U = 3 and W\ = 0.2 in Fig. 14. The 
behavior here is also quite interesting. The Lorenz num- 
ber is monotonic for total filling less than one, but then 
shows a large peak at a filling of 1. When the filling in- 
creases further, it develops a lower-temperature dip be- 
low a higher-temperature peak. This low-temperature 
dip can be advantageous to thermoelectric properties, 
but unfortunately the thermopower is too low (since it 
passes through zero) in this region to create a large ZT. 
Nevertheless, we do see interesting phenomena in the ZT 
curves as well. As we approach a total filling of 1, the 
low-temperature shoulder gets pushed closer to T = 0, 
and it remains at a fairly high value, above 1. When we 
move past 1, we find the low-T figure of merit is small 
due to the proximity to the sign change in S, but the 
high-temperature values are strongly enhanced. 

IV. CONCLUSIONS 

We have examined the thermal transport properties of 
the Falicov-Kimball model in a binary-alloy picture. By 
fixing the ion concentration as a function of temperature, 
we force the renormalizcd energy level of the localized 
particles to lie at the chemical potential as T — > 0, which 
is believed to be advantageous for thermal transport. In- 
deed, we find significant regions of parameter space with 
ZT > 1 at high temperature, and in the correlated in- 
sulators, we also find a small region of parameter space 
with ZT > 1 at low temperature. Our calculations ne- 
glect the lattice contribution to the thermal conductivity, 



which should have a limited effect on the high temper- 
ature calculations, but can destroy the low-temperature 
peaks in ZT if the lattice thermal conductivity is too big. 

We showed that generically, in a correlated insulator, 
one has S — > as T — > and C — > oo as T — > 0. Our 
analysis for the thermopower emphasized a relationship 
between S and the temperature dependence of the chem- 
ical potential, which appears to be general for systems 
that display a true gap. The generic dependences of S 
and L at low temperature imply that the generic ther- 
moelectric figure of merit would be small at low temper- 
ature. But in the infinite-dimensional hypercubic lattice, 
the noninteracting DOS is a Gaussian, which implies that 
the system really possesses a pseudogap in the correlated 
insulator, with the interacting DOS exponentially small 
in the gap region. This has a significant effect on the 
thermal transport because the relaxation time is not ex- 
ponentially small within the "gap region" , and in partic- 
ular, it can produce a low-temperature peak to ZT that 
moves to lower temperature as W\ approaches 0.5 and as 
U is tuned to lie close to the metal-insulator transition. 
If we took the hopping energy scale t* to be on the order 
of 1 eV, then the low-temperature peak in ZT can easily 
occur below 50 K [see Fig. 10 where the peak in panel 
(a) lies at about T = 0.003]. The key issue is whether or 
not the lattice thermal conductivity would wash out this 
effect. 

One question to ask, then, is can one find a way to in- 
troduce an exponentially small DOS into the gap region 
of a correlated insulator (if they do not appear in the 
bulk system)? The answer is yes, and it can be done by 
creating heterostructures of the correlated material and 
metals. The metallic DOS will "leak" into the correlated 
insulator, with a characteristic length scale, and create 
small subgap DOS within the system. Hence we believe 
the heterostructure idea of Mao and Bedell 7 or Rontani 
and Sham 6 can indeed allow one to get large peaks in 
the low-temperature figure of merit (if the exponentially 
small DOS also leads to a much larger "effective" scat- 
tering time). Since a heterostructure will also reduce 
the lattice thermal conductivity it is possible that the 
low-temperature peak could be realized, but it is likely 
to require a careful tuning of the correlation gap, and 
the thicknesses of the metallic and correlated layers of 
the heterostructure. One also has to be able to main- 
tain the correlated behavior within the thin layers of the 
heterostructure. 

Another potential way to create states within the gap 
is simply via a thermal rearrangement of the DOS within 
the bulk system as T increases. In the Falicov-Kimball 
model, the interacting DOS is T- independent, but in 
other correlated systems (like the Hubbard or periodic 
Anderson models) the correlated DOS does depend on 
T. It would be interesting to see if the creation of an ex- 
ponentially small DOS via thermal activation could allow 
for a peak in the low-temperature thermoelectric figure 
of merit (but this cannot be studied with the Falicov- 
Kimball model). 
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The situation at high temperature is much more 
promising. We find that generically both the electrical 
and thermal conductivities become flat at high tempera- 
ture, so the Lorenz number decreases like 1/T. The ther- 
mopower also becomes flat, with a magnitude growing as 
the system is made more particle-hole asymmetric. The 
net effect is a thermoelectric figure of merit that grows 
linearly in T, with a potentially large slope. We found 
the figure of merit usually improved at high temperature 
when the electronic filling was pushed higher in the band, 
and that there was no magical need to tune the electronic 
chemical potential to lie in the "gap region" , rather the 
enhancement was generic in correlated systems. 

We finally note that our scattering time t(u>) is never 
close in appearance to a delta function in this system. It 
can develop large asymmetry with a large peak lying at 
one side of the chemical potential, but the peak width is 
always broad, and determined by the effective bandwidth 
of the conduction electrons. This is, of course, because 
we have no hybridization in this model, which precludes 
the appearance of a sharp Abrikosov-Suhl resonance in 
the DOS and the similar formation of such a structure 
in t(u>). It would be interesting to see how the situation 
could change if hybridization was included, but this re- 
quires significantly more sophisticated numerical efforts. 
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